Inverse Freezing in Mean-Field Structural Glasses 
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A disordered spin model suitable for studying inverse freezing in fragile glass-forming systems is 
introduced. The model is a microscopic realization of the "random-first order" scenario in which 
the glass transition can be either continuous or discontinuous in thermodynamic sense. The phase 
diagram exhibits a first-order transition line between two fluid phases terminating at a critical point. 
When the interacting degrees of freedom are entropically favoured an inverse static glass transition 
and a double inverse dynamic freezing appear. 



Inverse melting and inverse freezing occur when a crys- 
talline or amorphous solid reversibly transforms into a 
liquid upon cooling. This unusual phase behaviour was 
first predicted by Tammann in 1903, and is generally 
considered to be rare because it does involve a counter- 
intuitive increase of thermal disorder as the temperature 
is loweredi. An example of biological relevance is pro- 
vided by elastirim. In the past few years "inverse tem- 
perature transitions" of this type have attracted a re- 
newed interest as they have been observed in a variety 
of soft matter systems including polymers^, colloids^ and 
micelles^. While the responsible physico-chemical inter- 
actions may depend on the system under study, it has 
been recognized that a large enough degeneracy of the de- 
grees of freedom interacting at low temperature provides 
a simple mechanism whereby inverse melting or inverse 
freezing may generally occur^. 

The idea is easily described by considering an ensemble 
of polymers that have a low temperature "folded" state in 
which they are mutually weakly interacting, and a higher 
temperature "unfolded" state which is favored entropi- 
cally and in which they interact strongly with each other. 
As temperature is increased, each polymer stretches out 
to reach the other polymers, the resulting entangling thus 
may lead to a glass transition. To obtain a minimal model 
of freezing by heating, one can consider^ spins taking val- 
ues 0, ±1, and a Hamiltonian consisting of a term of 
favouring the "folded" states = 0, and an interaction 
term ^j^- JijGiCfj that is active in the "unfolded" states, 
Ui = ±1. The entropic favouring of the latter is enhanced 
by making them r-fold degenerate. If the interactions 
matrix J is taken from the Gaussian ensemble, one then 
obtains a reentrant spin-glass phaseSiSiSiifliii. 

In this paper, the inverse freezing problem is addressed 
in the context of mean-field models of structural glasses 
by using the above mechanism of entropy-driven reen- 
trance. Several reasons make such a problem interesting. 
The "random first-order" scenario for the glass transi- 
tion^^ predicts that upon cooling fragile glass-forming 
liquids undergo a purely dynamic arrest before a ther- 
modynamic singularity occurs at a lower temperature 
(or higher density). The dynamic arrest is the relevant 
one from an experimental point of view and, in order to 
compare observations with theoretical predictions, one 
should consider the effect of degeneracy on the dynam- 



ics. The point is important because a reentrant glass 
transition has been recently predicted by mode-coupling 
theory^"^ and found in colloids^ and micelles with attrac- 
tive interaction^. The second and more general question 
that arises concerns the interplay of the two (static and 
dynamic) glass transitions and its effect on glassy be- 
haviour— when a reentrance in the phase diagram takes 
place. 

In order to answer the above questions, we consider a 
disordered system of N spin-1 variables with Hamiltonian 



H = —2 ^ JijffiUj 



= 0, ±1 (1) 



where J is a symmetric random orthogonal matrix (with 
Jii — 0); and D is a crystal field playing a role similar to 
the chemical potential: Increasing D will favour (7^ = 
states and reduce the effect of frustration. The case with 
binary spin variables (di = ±1) corresponds to the stan- 
dard Random Orthogonal Model (ROM) studied by Mari- 
nari, Parisi and Ritort — , which is known to be glassy 
at low temperatureiSiiSiii. It should be emphasized that 
quenched disorder is not crucial as ROM shares the same 
basic phenomenology with systems having deterministic 
interactionsASJ^ . 

The free-energy of the model (Q) can be evaluated as 
in the standard ROM by using the replica method and 
the identitjii^: exp(Tr JA) = exp(iVTr G(A/iV)), where 
A is a symmetric matrix of finite rank, the overbar is the 
average over the quenched disorder which is defined by 
the Haar measure on the orthogonal group, and 



G(z) = ^ 2 In . (2) 

Averaging the replicated partition function gives 

/oo 
l[dAabdQabexp{-/3f[Q,A]) , (3) 

a,b 

where a, 6 = 1, . . . , n are replica indexes, and the Parisi 
order parameter Qab = {(JaCb) includes the diagonal 
terms, which corresponds to the density of ±1 spins, 
Qaa — P- The free energy f[Q,A] reads 



-/3/[Q,A] = iTrG(4/3Q) 



Tr(AQ) + lnZo, (4) 
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FIG. 1: (Color online) Phase diagram of ROM spin-1 in the 
temperature-crystal field plane. The continuous (red) line 
is the dynamic freezing, while the dashed (green) line is the 
static glass transition. The dotted (blue) line is the first-order 
transition between the two paramagnetic (PMi and PM2) 
phases. Spinodals are shown as dotted (violet) light lines. 
The full dot is the critical point terminating the coexistence 
line. The diamond symbol is the tricritical point separating 
the continuous and discontinuous glass transition (at lower 
temperature). 

where is the single-site partition function 

Zo = ^ CXp ^ Kab TaCTb ^ PDJ2<^1\ ■ (5) 

{o"a} \a.,b a J 

In order to proceed one now needs to specify an ansatz for 
Q and A, and then consider the zero-replica limit n — > 0. 
The fluid-fluid transition. The simplest case is the 
replica symmetric ansatz: Qab ~ {p ~ q)Sab + Q, ^ab — 
(/i — X)8ab + A, where b is the Kronecker symbol. In this 
approximation the free energy reads: 

- jDz\n(l + 2e''-^-'^^ cosh(z\/2A)) , (6) 

where Dz = dze~^ ^^/^/2n and q, p, A and p., are self- 
consistently determined by the saddle-point equations. 
At sufficiently large temperature/crystal field, q = X = 0, 
and one recovers the annealed free energy 

/3/a„„ = -^G{Af3p)+f3Dp-soip), (7) 

where so{p) = —(1 — p)log(l — p) — plogp + plog2 is 
the entropy of a noninteracting spin-1 system, and the 
density p = p{P, D) satisfies the implicit equation 

f3D = log ^iizf) + 2/3G'(4/3p) . (8) 
P 

Equation (jHJ exhibits multiple solutions. There is an un- 
stable phase (with negative susceptibility) and two para- 
magnetic fluid phases (PMi and PM2) between which a 



FIG. 2: (Color online) Energy density e vs crystal field D 
at fixed temperature T, in a slow crunching MC experiment. 
The lines are the analytic results of the annealed approxima- 
tion. 

first-order transition occurs. The latter terminates at 
a critical point located at — 2.644, ~ 0.5506, 
see Fig. ^ The critical point does not appear in spin- 
1 REM (Ref. Ill), showing that the two models are not 
equivalent in the large- Z? regime. Along the first-order 
transition line two distinct phases having equal free en- 
ergy coexist, and in both it is possible to go continuously 
around the critical point from one coexisting phase to the 
other by appropriately varying D and T. This first-order 
transition is a general feature of spin-1 and lattice-gas 
systems with disordered (Gaussian or orthogonal) inter- 
actions, though it is sometimes missed. For D > 2 the 
density of ±1 spins decreases upon cooling. This be- 
haviour is somehow unexpected if compared to what hap- 
pens in the Ghatak-Sherrington model, i.e. Eq. Q with 
Gaussian disorder. In the latter, the density increases 
upon cooling and the critical point is absent. 
The static and dynamic glass transition. At low 
temperature/crystal field, the annealed entropy, 

Sa„„ = soip) + ^G{Af3p)-2PpG'{APp), (9) 

becomes negative, suggesting that replica symmetry has 
to be broken. Within the one-step replica-symmetry- 
breaking ansata2£, Q and A are block diagonal matrices, 
where the blocks have size m x m. Inside the blocks 
Qab = [p- q)5ab + q, Aah = (/i - X)5ab + A. Then the free 
energy becomes 

1 — 77? 

= Gmp~q)) + \q{m~l) 

-^G{AP{p-q + qm))+ tip (10) 
2m 

-^In Jdz (l + 2e^-^-'^^cosh(zV2A)y" . 

Expanding near m = 1 gives /3/msB — /3/ann ~ {fn — i)V, 
that allows to locate the static and dynamic transi- 
tion through the effective potential V — — /3 '^■^^"^^ 
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FIG. 3: (Color online) Energy density e vs. temperature T at 
fixed crystal field D, in a slow cooling MC experiment. The 
lines are the analytic results of the annealed approximation. 



fRefs. IT^ and IT^ . The two glass transition lines Ts{D) 
and Td(-D) are shown in Fig. ^ In the Ising spin limit, 
D ~oo, we get Tg ~ 0.26 and — 0.535, consis- 
tently with Refs. ITsI and IT^ Several interesting features 
can be observed, (i) The temperature at which the an- 
nealed entropy vanishes is very close to Ts{D) for any 
value of D, meaning that the glassy phase of our model 
is similar to that of spin-1 REM (Ref. Hil). (ii) There is 
a tricritical point at T* ~ 0.036, D* ~ 2.0 below which 
the nature of the glass transition changes from second 
to first order in the thermodynamic sense. No apprecia- 
ble irreversibility effects are observed across the second- 
order glass-transition line when D loops start from the 
PMi phase. Fig. |2] (main frame) , whereas there is latent 
heat and hysteresis across the first-order glass-transition 
line, see Fig. |21 (inset), (iii) The dynamic freezing line 
penetrates the PM2 phase up to the spinodal line (at 
D ~ 2.1), see Fig. For D < 2.1 the system is dy- 
namically unable to reach equilibrium at low tempera- 
ture even in the case in which the ground is trivial, that 
is for 2 < D < 2.1, see Fig.|31 (iv) Crunching the system 
from high to low D at fixed temperature, Td > ^ > Tg, 
leads to a purely dynamic freezing with no underlying en- 
tropy crises. Monte Carlo (MC) results shown in Figs. [21 
and 131 are for a system of size TV = 512 and a very slow 
annealing rate (lO*" MC sweeps per unit variation of D 
and T, respectively). 

Inverse temperature glass transitions. To take 
into account inverse freezing phenomena, the interact- 
ing states, (Ti = ±1, are now given an entropic advan- 
tage by a degeneracy ratio r > 1 with respect to the 
noninteracting states, Ci = 0. The opposite case in 
which r < 1 will not be discussed here: That would 
lead only to a reentrancc of PMi within the PM2 fluid 
phase, which, for our purpose, is less interesting. One 
can easily see that including degeneracy in the system 
corresponds to changing the crystal field in Eqs. (6-10) 
as follows: D ^ D — Tlogr (Ref. |^. In fact, the depen- 
dence of the effective potential on the degeneracy ratio 
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FIG. 4: (Color online) Phase diagram of ROM spin-1 for de- 
generacy ratio r = 50. Spinodal lines are not shown here. 



enters only through the density variable, p, which satis- 
fies Eq. (|H1). Having the meaning of the different phases 
already clarified, we need only to investigate how the 
phase boundaries in Fig.^are modified accordingly. The 
analysis of the effective potential shows that for small r 
the reentrancy effect is weak, while for large enough val- 
ues of r a rather interesting reentrant behaviour appears. 
We discuss the results for the case r = 50, which corre- 
sponds to the phase diagram presented in Fig. ^ First 
of all, we see that the two glass-transition lines never 
cross each other and that upon heating the inverse static 
glass transition occurring in the range 2 < D < 2.18 is 
anticipated by inverse freezing. This is due to the pen- 
etration of the dynamical arrest line in the PM2 fluid 
phase. At intermediate crystal field, 2.18 < I? < 2.37, 
there is inverse freezing without a static glass transition. 
This prevents the application of the Kauzmann paradox, 
since in this range of crystal field the equilibrium phase 
at low temperature is the PM2 fluid. Finally, at larger 
crystal fleld, 2.37 < D < 2.58, there is a double inverse 
dynamical freezing (again with no underlying static glass 
transition) that goes through the PMi and PM2 fluid 
phases (i.e., on cooling one would observe the sequence 
of transitions PM1-G-PM1-G-PM2). In this region of the 
phase diagram the packing density of ± 1 spins decreases 
with the temperature, i.e., the glass state at lower tem- 
perature is less dense than the one at higher tempera- 
ture. This behaviour is reminiscent of that predicted by 
mode-coupling theorjii^ and observed in experiments on 
attractive colloids^ and micelles^. 

Conclusions. To summarize, we introduced a gen- 
eralisation of ROM (Ref. llEI) allowing investigation of 
inverse freezing phenomena in fragile glass-forming liq- 
uids through a mechanism of entropy-driven phase reen- 
trancc®. The model is a microscopic realization of the 
"random first-order" scenario for the the structural glass 
transitions^. However, in our case the glass transition 
can be either continuous or discontinuous in thermo- 
dynamic sense. This is due to the presence of a first- 
order transition between two fiuid phases. Similar re- 
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suits are obtained by using purely biquadratic interac- 
tions, Jijafa"^, that would correspond to a lattice-gas 
ROM {af — > Hi = 0,1). Notice that our results have 
no counterpart in the quantum version of ROM studied 
in Ref. Rather, they have some resemblance with 
those obtained in Refs. and reproducing a num- 
ber of experimental observations on the dipolar spin-glass 
LiHoa;Yi_2:F4 in external fieldS^. In fact, the present 
model can be considered as the insulating limit of a gen- 
eralisation of the itinerant electron model (see Ref. 12711 
with random orthogonal interactions. 

At large enough degeneracy an inverse static glass tran- 
sition and a double inverse dynamic freezing occur. The 
latter reproduces qualitatively some features observed in 
recent experiments on colloidal and copolymer-micellar 
systems with short-range attractioniiS,. Nevertheless, the 
possibility of describing the glass-glass transitioniiSiiS^ 



in the present setting remains unclear (see, however, 
Rels. andl29ll and deserves further investigation. Let 
us finally mention that this work can be extended by 
considering both bilinear and biquadratic interactions, 
similarly to what has been done in Ref. 10. One can 
also include a three-body interaction term to mimick mi- 
croemulsion (see, e.g., Refs.'s^). Depending on the rela- 
tive strength of these interactions and the nature of the 
quenched (Gaussian versus orthogonal) disorder an even 
richer variety of phases is expected. That could lead to a 
better understanding of the glassy behaviour of complex 
liquids and soft matter systems. 



It is a pleasure to thank L. Cugliandolo, D. Dean, and 
especially J. Kurchan for discussions and suggestions. 
Support of the EVERGROW project is acknowledged. 



^ F.H. Stillinger and P.G. Debenedetti, Biophys. Chem. 105, 
211 (2003); F.H. Stillinger and RG. Debenedetti and T.M. 
Truskett, J. Rhys. Chem. B 105, 11809 (2001) 

^ D.W. Urry, J. Rhys. Chem. B 101, 11007 (1997); L. De- 
belle and A.M. Tamburro, Int. J. Biochem. Cell Biol. 31, 
261 (1999) 

^ C. Chevillard and M. Axelos, Colloid Rolym. Sci. 275, 537 
(1997); S. Rastogi, G.W.H. Hohne and A. Keller, Macro- 
molecules 32, 8897 (1999) 

K.N. Rham, A.M. Puertas, J. Bergenholtz, S.U. Egelhaaf, 
A. Moussaid, R. N. Pusey, A.B. Schofield, M.E. Gates, M. 
Fuchs, and W.C. Poon, Science 296, 104 (2002); T. Eckert 
and E. Bartsch, Rhys. Rev. Lett. 89, 125701 (2002) 
^ S.H. Chen, W.R. Chen and F. Mallamace, Science 300, 
619 (2003) 

^ N. Schupper and N. M. Shnerb, Rhys. Rev. Lett. 93, 
037202 (2004); Rhys. Rev. E 72, 046107 (2005) 
A. Crisanti and L. Leuzzi, Rhys. Rev. Lett. 95, 087201 
(2005); Rhys. Rev. B 70, 014409 (2004) 

* S.K. Ghatak and D. Sherrington, J. Phys. C 10, 3149 
(1977); P. Mottishaw and D. Sherrington, J. Rhys. C 18 
5201 (1985) 

^ F.A. da Costa, C.S.O. Yokoi and S.R. Salinas, J. Phys. A 

27, 3365 (1994); F.A. da Costa, F.D. Nobre and C.S.O. 

Yokoi, J. Phys. A 30, 2317 (1997) 
^° J.J. Arenzon, M. Nicodemi and M. Sellitto, J. Phys. I 6, 

1143 (1996); M. Sellitto, M. Nicodemi and J.J. Arenzon, 

J. Phys. I 7, 945 (1997) 
" G.R. Schreiber, Eur. Rhys. J. B 9, 479 (1999) 
^2 T. R. Kirkpatrick and P. G. Wolynes, Phys. Rev. B 36, 

8552 (1987); T. R. Kirkpatrick, D. Thirumalai and P. G. 

Wolynes, Phys. Rev. A 40, 1045 (1989) 

K.A. Dawson, G. Foffi, M. Fuchs, W. Gotze, F. Sciortino, 

M. Sperl, P. Tartaglia, Th. Voigtmann, and E. Zaccarelli, 

Phys. Rev. E 63, 011401 (2000) 

L.F. Cugliandolo, in Slow relaxation and nonequilibrium 



dynamics in condensed matter, J.-L. Barrat et al. eds. 
(Springer, Berlin, 2003) 

E. Marinari, G. Rarisi, and F. Ritort, J. Rhys. A 27, 7647 
(1994); J. Phys. A 27, 7615 (1994) 

F. Rao, A. Crisanti and F. Ritort, Europhys. Lett. 62 869 
(2003) 

R. Cherrier, D.S. Dean and A. Lefevre, Phys. Rev. E 67, 
046112 (2003) 

J.R Bouchaud and M. Mezard, J. Rhys. I 4, 1109 (1994); 

W. Krauth and M. Mezard, Z. Rhysik B 97 127 (1995) 

P. Mottishaw, Europhys. Lett. 1, 409 (1986) 

M. Mezard, G. Parisi and M. Virasoro, Spin Glass Theory 

and Beyond, World Scientific (Singapore, 1987). 

A.M. Puertas, M. Fuchs and M.E. Gates, Rhys. Rev. Lett. 

88, 098301 (2002) 

F. Sciortino, P. Tartaglia and E. Zaccarelli, Phys. Rev. 
Lett. 91, 268301 (2003) 

F. Ritort, Phys. Rev. B 55, 14096 (1997) 

Y.Y. Goldschmidt, Rhys. Rev. B 41, R4858 (1990); Th.M. 
Nieuwenhuizen and F. Ritort, Rhysica A 250, 8 (1998) 
L.F. Cugliandolo, D.R. Grempel and C.A. da Silva Santos, 
Phys. Rev. Lett. 85, 2589 (2000); Phys. Rev. B 64, 014403 
(2001) 

W. Wu, D. Bitko, T.F. Rosenbaum, and G. Aeppli, Phys. 
Rev. Lett. 71, 1919 (1993) 

H. Feldmann and R. Oppermann, J. Rhys A 33, 1325 
(2000); LP. Castillo and D. Sherrington, Rhys. Rev. B 72, 
104427 (2005) 

A. Caiazzo, A. Coniglio and M. Nicodemi, Phys. Rev. Lett. 
93, 215701 (2004) 

A. Crisanti and L. Leuzzi, Phys. Rev. Lett. 93, 217203 
(2004); Phys. Rev. B 73, 014412 (2006) 

G. Gompper and M. Schick, Rhys. Rev. Lett. 62, 1647 
(1989); Phys. Rev. B 41, 9148 (1990) 



